Cooling rates were calculated for neutron stars of about one solar mass and 10 km radius, with magnetic fields from zero to about 10 to the 14th power gauss, for extreme cases of maximum and zero superfluidity. The results show that most pulsars are so cold that thermal ionization of surface atoms would be negligible. Nucleon superfluidity and crystallization of heavy nuclei were treated quantitatively, and more realistic hadron star models were chosen. Cooling rates were calculated for a stable hyperon star near the maximum mass limit, a medium weight neutron star, and a light neutron star with neutron-rich heavy nuclei near the minimum mass limit. Results show that cooling rates are a sensitive function of density. The Crab and Vela pulsars are considered, as well as cooling of a massive white dwarf star.